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THE  BEARING  CAPACITY  OF  FLOATING  ICE  PLATES 
SUBJECTED  TO  STATIC  OR  QUASI-STATIC  LOADS 
A Critical  Survey 

by 

Arnold  D.  Kerr 


INTRODUCTION 

Frozen  lakes  and  rivers  have  been  utilized  since  early  times  for  transportation  and  storage 
purposes.  In  Russia,11’  in  the  absence  of  bridges,  railroad  tracks  have  been  placed  over  frozen 
rivers  since  about  1890.  Floating  ice  plates  have  been  increasingly  utilized  as  airfields  for  the 
landing  of  aircraft,1  “ 18  1,8  ,,,  as  platforms  for  storage  in  logging  operations,1'  111  as  platforms 
for  the  construction  of  river  structures.*4  141  as  off-shore  drilling  platforms  in  the  northern  regions.10 
and  as  aids  in  various  other  civilian  and  military  operations.*  The  successful  defense  of  Lenin- 
grad during  World  War  II  was  greatly  facilitated  by  the  "ice  road"  ove.  Lake  Ladoga."  The  recent 
oil  discoveries  in  northern  Alaska  have  increased  the  interest  in  the  arctic  ice  cover  Tor  off-shore 
drilling  purposes.  A rational  utilization  of  floating  ice  plates  for  all  these  activities  requires  the 
knowledge  of  their  bearing  capacity  when  they  are  subjected  to  loads  of  short  and  long  duration. 
Such  information  is  also  needed  for  the  design  of  icebreakers.48  *" 

Field  observations  reveal  that  when  a vehicle  is  small  and  rel;  tively  heavy  it  may  break 
through  the  ice  plate  immediately  after  placement.  In  such  cases,  he  plate  response  may  be  con- 
sidered elastic  up  until  failure.  For  relatively  light  vehicles,  the  ice  plate  deforms  elastically 
at  the  instant  of  loading,  but  sustains  the  load.  However,  as  time  progresses,  the  ice  plate  con- 
tinues to  deform  in  creep,  especially  in  the  vicinity  of  the  vehicle,  and  after  a certain  time  inter- 
val the  vehicle  may  break  through  the  ice. 

In  the  past,  numerous  attempts  have  been  made  to  determine  the  bearing  capacity  of  floating 
ice  plates  subjected  to  vertical  loads.  Particularly,  since  World  War  11,  many  papers  containing 
test  data  and  related  analyses  have  been  published.  However,  in  spite  of  these  publications, 
there  is  as  yet  no  reliable  analytical  method  for  predicting  the  bearing  capacity  of  floating  ice 
plates  subjected  to  static  or  dynamic  loads.  This  is  particularly  the  case  for  floating  ice  plates 
reinforced  by  pressure  ridges,  a phenomenon  often  encountered  in  the  Arctic,’1  ,,B  for  which  not 
even  test  data  can  be  located  in  the  literature. 

One  of  the  main  reasons  for  the  lack  of  reliable  methods  for  determining  the  breakthrough 
loads  of  ice  plates  is  the  difficulties  introduced  by  the  fact  that  the  lower  surface  of  an  ice 
plate  is  always  subjected  to  the  melting  temperature  of  about  0°C,  at  which  the  mechanical  prop- 
erties of  ice  vary  drastically  with  small  changes  of  temperature.  Other  difficulties  are  the  depend- 
ence of  tbs  mechanical  properties  of  the  ice  plates  upon  the  rate  of  freezing,  the  velocity  of  the 


* Refs.  5,  1 7 , 25,  40,  82,  120. 
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Uiilci  below  th<-  plate  during  the  fiee/ing  pioeess.  the  salmitv  of  the  walei,  eft  Diseussums  of 
tin  mechanical  properties  of  we  have  terentlv  heen  presented  by  Voit kovsk i i . ,4'  Weeks  and 
Assnr,1*’  ,4'  Lavrov,”  and  llugorodskn  el  at.1’ 

Another  main  reason  is  the  lai  k ol  effective  oomniunieal ion  among  the  various  investigators, 
p.ntly  eaused  by  the  landmine  hamei.  This  has  iesulled  m the  duplieation  of  analyses  and  tests, 
olten  tendeied  useless  by  the  same  shortcomings.  Also,  the  intiodurt ion  of  1111  oneet  solutions  foi 
Moating  iee  plates  and  their  suhsei|iient  utilization  for  eomparison  with  test  data  have  not  helped 
in  solving  the  pi  oh  le  ins  under  eonsideialion. 

The  purpose  ol  this  report  is  to  present  a eritieal  survey  of  the  literature  on  the  hearing 
capa<  itv  of  floating  iee  plates.  Tnst,  the  various  analytical  attempts  to  determine  the  bearing  ca- 
Ptieitv  are  leviewed,  grouped  aeeoiding  to  the  used  ‘'failure  eriteuon."  Tli.s  is  followed  hv  a dis- 
cussion ol  test  data  and  then  i elation  to  the  analytical  results.  The  tepoit  coni  ludes  with  a sys- 
teniaiu  summary  of  insults,  a discussion  of  otiseived  short  coinings,  and  teeommendat  ions  lor  needed 
investigations.  It  is  hoped  t!  it  this  smvey  and  summary  of  results  will  establish  a sense  of  direc- 
tion in  the  investigations  and  will  lonliibute  lowaid  developing  methods  for  deteiinining  the  hearing 
• apat  itv  ol  floating  iee  plates. 


ANALOGY  METHOD 


The  analogy  method  of  predicting  the  bearing  capacity  of  a floating  ice  plate  subjected  to 
a sialic  vertical  load,  discussed  by  Kornilov,”  is  based  on  the  notion  of  the  analogy  of  two 
plates.  Koiunov  assumed  that  tin*  ice  plates  under  consideration  are  homogeneous  and  isotropic 
and  that  for  two  plates  with  thicknesses  ft,  and  h2  the  corresponding  failure  moments  M,  and  U2 
in  cylindrical  bending  are 


U 


t 


u 


(1) 


Assuming  that  the  failure  stress  <;f  for  the  two  plates  is  the  same,  it  follows  that 


M,  ft.* 


(2) 


Considering  the  eflecl  of  two  different  loads,  P(  acting  on  the  plaie  with  thickness  ft,  and  P2  act- 
ing on  the  plate  w it h thickness  h2,  Korunov  assumed  that  U is  proportional  to  P,  and  obtained 
from  eq  2 


P, 


Cl) 


Equation  3 may  he  rewiitten  as  follows 


Note  lliat  ’ i)  •'!  was  usect,  in  I9'I8,  l>y  Mnskalov  Oct.  HC , p.  SI). 


Tilt,  HEARING  CAPACITY  Oh'  FLOATING  ICE  PLATES 


all 


Ah 


3 


(1) 


wlK're  ^ * P>  According  to  tin*  above  met  lux),  if  an  allowable  load  P , of  an  ice  plate  of  thick* 

ness  h , is  known  tfio.ii  a test),  then  the  allowable  load  P „ of  an  tee  plate' of  different  thic  kness 
ntav  be  computed  tt  the  „f  values  are  the  same  for  both  plates,  rims,  the  coefficient  A in  eq  1 is 
to  lie  determined  liom  a s|>et  if  it  test. 


Souu  shortcomings  in  the  derivation  of  eq  1 weie  discussed  l>v  Lagutin  and  Slniliua».Tft  It 
should  also  lie  noted  that  in  a Hunting  ice  plate  the  bending  stress  distribution  may  not  he  linear 
across  the  plate  thickness  therefore,  eq  1 in  the  above  derivations  may  not  lx*  admissible. 
Nevertheless,  because  of  its  extreme  simplic  ity  and  its  agreement  with  various  test  results,  eq  1 
found  wide  popularity,  as  shown  in  the  followup:  table  (valid  for  P ut  metric  tons  and  h in 


centimeters) 


Source Load A 

Kotuiiov  ' 

lVschausku"'  OBI 

Letx-dev" 

Zubov1"  0.0166 

Instructions  of  the  Kip;  nicer  me  Whet  led  vehicles  0.0070 

Committee  of  the  Red  Army"  Tracked  vehicles  0 0123 

1 ysukhin"  Wheeled  vehicles  0.0082 

Tracked  vehicles 0.0123 


To  demonstrate  the  use  of  eq  1 let  us  deteimine  the  necessary  ice  thickness  for  the  crossing 
of  a river  by  a truck  weighing  30  metric  tons,  accoidmg  to  Korunov."  Using  eq  1 the  necessary 
ice  thickness  is 

h \ TOO  ^ Tq  iq  , f,  (,o  t.|u. 

Additional  examples  of  the  use  of  eq  1 were  presented  by  Moskatov.*"  l.ysiiklun* 2 and  Gusev.*' 

In  order  to  take  into  consideration  the  olfects  of  temperature,  the  dimensions  of  load  distri- 
bution. and  the  salinity  of  ice,  Zubov1’4  modified  eq  1 as  follows. 

Palt  KM  SA  h~  (f,) 

where  K,  M.  and  s are  the  corresponding  correction  coefficients.  Discussions  of  this  extension 
are  presented  in  ref.  7f>  and  151. 

B lsing  Ins  woik  on  field  experience  with  fresh  water  ice,  Korniiov,70  in  1956,  modified  eq  4 
by  introducing  a correction  coefficient  n which  takes  into  consideration  the  condition  of  the  ice 
as  follows 


P 


att 


n 


in  tons. 


(6) 


In  the  above  formula  A 


0.01  and  n is  related  to  <>(  as  follows 
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10 


20  30 

a.  (kg/cm2) 


i 

40 


Figure  1.  Correction  coefficient  n as  related 
to  failure  st.  s of  (Korunov10). 


Of  (kg/cm*) 

5 

12 

17 

25 

38 

n 

4.8 

2.0 

1.4 

1.0 

0.6 

for  T<  -7  C. 


A graph  of  these  v titles  is  given  in  Figure  1.  This 
graph  may  be  represented  by  the  equation 


25 

n = — . 


Substituting  this  into  eq  6,  we  obtain  for  T < -7  C 


P„,,  «•  o.h ^ in  tons 

a11  2500  r 


or 


P II  = 0,4af/i2  in  kilograms.  (7) 


ot  values  were  stipulated  by  Korunov”  for  five  types  of  ice.  Korunov’0  also  introduced  another 
correction  coefficient  for  thaw  temperatures. 


Ml  TH0D  BASED  ON  THE  BENDING  THEORY  OF  ELASTIC  PLATES 
AND  THE  CRITERION  e of 


This  me. hod  of  predicting  the  hearing  capacity  of  a floating  plate  subjected  to  loads  of  short 
duration  consists  of  the  following  three  steps: 

1.  Determination  of  the  maximum  stress  amax  in  t he  floating  ice  plate  due  to  a given  load,  as- 
suming that  the  ice  plate  is  elastic. 

2.  Determination  of  the  load  at  which  the  first  crack  occurs  Pcr,  utilizing  the  critirion 


”max  °t  • 


(8) 


3.  Correlation  of  Pcr  with  the  breakthrough  load  Pf.  This  step,  disregarded  by  many  investi- 
gators, is  needed  because,  according  to  field  tests,  for  various  plate  geometries,  the  occurrence 
of  the  first  crack  does  not  cause  breakthrough;  therefore,  for  these  cases  Pf  > Pcr. 

In  the  criterion  in  eq  8,  «f  is  the  “failure  stress,"  It  is  usually  obtained  by  loading  a floating 
ice  beam  to  failure  and  then  computing  the  largest  bending  stress  at  which  it  failed.  In  the  located 
literature,  is  determined  using  the  classical  bending  theory  of  thin  elastic  plates.  Tliese  results 
are  reviewed  in  the  following. 

The  response  of  a homogeneous  and  Isotropic  elastic  plate  that  rests  on  a liquid  and  is  sub- 
jected to  a static  vertical  load  q is  described  by  the  partial  differential  equation 


where 


D w ♦ yw  q 


w(x,y)  = plate  deflection  at  (x,  y) 

D ■=  flexural  rigidity  of  the  plate 
y = specific  weigh*  of  the  liquid. 


(9) 
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Figure  2.  A flouting  ice  plate  subjected  to  a distributed 
loud  q over  a circular  area  of  radius  a. 


0 60' 


0 50'  \ 


040' 


C a) 


0 30 


0 20' 


Solutions  for  the  infinite  plate  subjected  to  a con- 
centrated load  P,  and  to  a load  uniformly  distributed  over 
a circular  area,  were  presented  by  Hertz*1  in  1881.  In 
1929  Bernshtein*  utilized  this  discussion  for  the  deter- 
mination of  the  allowable  load  for  an  infinite  ice  plate. 
Using  the  criterion  in  eq  8,  in  conjunction  with  the  solu- 
tion for  an  infinite  plate  subjected  to  a uniform  load  over 
a circular  area,  as  shown  in  Figure  2,  Bernshtein  obtained 


o,  h* 


cr  1lUe)C(n)  r 


(10a) 


o 10- 


where 


O 02  * 04  ’ 36  08 

a 

Figure  3.  C(a>  vs  a graph.* 


i-  = Poisson's  ratio  for  the  plate  material 
C(a)  a given  function  of  a a //,  as  shown  in  Fig- 
ure 3 

a rudius  of  the  circular  area  subjected  to  the 
uniform  load  q P/rra 2 

i <577 
o E*Vtiai 


If  omax  = Of  is  a valid  criterion,  then  Pcr  is  the  load  intensity  at  which  the  plate  cracks. 

To  demonstrate  the  use  of  eq  10a,  Bernshtein  computed  the  omaj;  due  to  a railroad  car  weighing 
24  tons  for  a 70-cm-thick  ice  plate  as  follows  (ref.  8,  para.  20) 

Assuming  that  E = 550,000  t/m*  and  v = he  obtained 


11.50  m. 


He  then  assumed  that  the  effect  of  the  weight  of  the  railroad  car  may  be  represented  by  a load  uni- 
formly distributed  over  a circular  area  with  radius  a 1.54  m.  Hence,  a - a /L  = 0.134.  From  Fig- 
ure 3,  it  follows  that  C(a)  0.417.  For  the  above  values  eq  10a  yields 


amax 


24,000  x 3x-x0.417 
1 

(70) 2 


8. 16  kg/cni2. 


* This  is  a modified  graph.  In  the  original  version,*  the  C(a)  presented  is  for  P in  Ions,  h in  meters,  and  a 
in  kilograms  per  square  centimeter. 
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The  next  step  is  to  check  whether  <’|nax  Additional  numerical  analyses  are  Riven  in  ref  8. 

Other  numerical  examples,  based  on  the  Berushteiu  solution,  were  presented  by  Volkov1"  in  1940 
and  by  Bregman  and  Prosknriakov  (ref.  11,  part  IV,  section  7)  in  1943. 

The  determination  of  the  load  Pfr  for  a flouting  infinite  p lute  based  on  eq  9,  the  criterion  in 
eq  8,  and  the  assumption  that  the  load  q P/(nu2)  is  distributed  uniformly  over  a circular  region 
of  radius  a was  also  presented  by  Wyman1”  in  1930,  Kubo”  in  1958,  and  Savel'ev”1  in  1963.  Wyman 
obtained  for  the  load  Pcr  the  equation 


v a 

3(1  * i ) kei ' a 


(10b) 


This  is  identical  to  eq  10a,  noting  that 


C(a) 


kei  '(q) 
rra 


(ID 


The  determination  of  Pc[,  assuming  that  the  uniform  load  is  distributed  over  a square  area  with 
sides  b,  was  obtained  by  Golushkevich”  in  1944.  The  derived  expression  yields  loads  that  are 
very  close  to  those  obtained  from  eq  10. 

Solutions  for  an  infinite  plate  were  also  presented  by  Schleicher1"  in  1926,  Korenev*0  in  1954, 
Korenev*1  in  I960,  and  Korenev  and  Chernigovskaia”  in  1962. 

A solution  for  the  infinite  plate  subjected  to  a row  of  equidistant  loads  was  piesented  by 
Westorgaard1”  in  1923,  in  terms  of  a trigonometric  series.  Solutions  to  similar  problems  (periodic 
load  distribution),  also  in  terms  of  trigonometric  series,  were  presented  by  Lewe*0  in  1923,  Mhller”  ** 
in  1952,  and  Panfilov100  10’  in  1903  and  1964.  Sh -kilter  and  Vinokurova1"  discussed  related  prob- 
lems in  1936. 

Since  eq  9 is  linear,  it  appears  that  when  the  plate  is  subjected  to  several  loads,  the  method 
of  superposition  should  be  used.  This  idea  was  demonstrated  by  Kerr50  in  1959  for  the  solution  of 
the  floating  ice  plate  subjected  to  a row  of  equidistant  loads.  A major  advantage  of  this  approach 
is  that  the  distribution  of  the  loads  on  the  floating  plate  may  be  arbitrary,  whereas  the  use  of 
trigonometric  series  is  suitable  only  when  the  loads  act  along  straight  lines,  all  loads  along  a line 
are  of  the  same  intensity  and  distribution  and  the  distance  between  them  is  the  same. 

The  analysis  of  floating  ice  plates  for  arbitrary  load  distributions  may  be  greatly  simplified 
by  utilizing  influence  surfaces. 1,1  Charts  of  such  surfaces  were  presented  by  Pickett  and  Ray11* 
in  1951  for  concrete  pavements.  Influence  surfaces  for  bending  moments,  more  suitable  for  ice 
plate  problems,  were  presented  by  Palmer’*  in  1971.  Palmer’s  charts  could  also  be  used  for  the 
determination  of  load  distributions  on  the  plate  that  yield  the  largest  possible  tending  moments. 

An  attempt  to  solve  such  a problem  without  influence  surfaces  was  made  in  1965  by  Nevel  and 
Assur.’*  They  considered  the  problem  of  the  most  unfavorable  distribution  of  crowds  on  a floating 
ice  plate  from  the  point  of  view  of  bearing  capacity,  based  on  the  criterion  in  eq  8.  This  problem 
was  recently  analyzed  by  Palmer"*  using  influence  surfaces. 

Bernshtein’s*  eq  10a  is  shown  as  the  solid  line  in  Figure  4.  Shulman111  in  1946  simplified 
eq  10a  by  replacing  the  curve  for  0.07  < « < 0.65  with  a straight  line  described  by  the  expression 

Prr  „ 0.375  af(hZ  i 7.8a^  h *) 


(12a) 
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Figure  4.  Vc/<>{  lr  vs 


Panfilov'’  in  1960  proposed  the  expression 
P,f  0.376(1  , 1. 1 a )tt(h^ 


(12b) 


based  on  the  idea  of  a straight-line  approximation.  Panfilov's  approximation  (eq  12l>)  is  the  same 
as  the  one  presented  by  Shulnian,"1  since  for  v 0.3 


I ^ i?  '‘F. 


(13) 


Panfilov”  also  proposed  the  following  approximation 


Pcr 


3(1  4-  v)(0.6H2  f 0.019a2  - In  a) 


■a{hL. 


(11) 


However,  since  this  is  not  much  simpler  than  the  exact  expression  (eq  10a  or  10b),  its  usefulness 
is  questionable. 

In  1961,  Panfilov1”  attempted  to  derive  another  approximate  expression  for  Pcr,  assuming  that 
the  deflections  of  a floating  ice  plate  subjected  to  a concentrated  force  P may  be  expressed  ap- 
proximately as  follows: 


w(x,  y)  = wQ  eXp|--^(x  +y)|^sin  t-  cos  ^?j|sin  ^ c cos  — j 


where 


A 

From  the  equilibrium  equation 


(15) 


(16) 
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OC  9C 
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dx  dy. 


(17) 
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Panfilov  determined  the  only  unknown,  w0,  as 


K, 


P 


(18) 


Comparing  the  resulting  w(x,  y)  with  the  exact  solution,  and  finding  t hat  the  agreement  was  relatively 
dose,  Panfilov  determined  the  bending  moments,  using  the  relations  (ref.  142,  p.  81) 


WJ x,  y) 


-bKZ  ■ . 

VrJx^  <)y  2 ) 

M ix,  y)  d(—  ♦ „ °±) 

• U2  ,lxs  I 


(19) 


and  the  approximate  tv(x,  y)  given  in  eq  15.  For  the  bending  moments  under  the  load  P.  he  obtained 

(I  * i)P 


MAO,  0)  MAO,  0)  = 

* > R 


(20) 


Equating  this  expression  with  MCf  = nfh 2 6,  Panfilov  obtained  for  e the  expression 


pc, 


(21) 


At  this  point,  note  that  the  relative  closeness  of  the  approximate  and  exact  deflections  (in  the 
sense  of  comparing  two  graphs)  does  not  imply  that  the  second  derivatives  are  also  close.  Thus, 
for  example,  whereas  ttie  exac  t solution  for  the  classical  plate  theory  (used  by  Panfilov)  yields 
infinite  moments  under  the  concentrated  load  P*  Panfilov’s  approximate  solution  yields  the  finite 
value  shown  in  eq  20.  This  point  may  be  demonstratec  further  by  comparing  the  graphs  for  the 
bending  moment  Mx(x,  0)  based  on  eq  15  and  or  the  exact  solution.  It  may  be  shown  that,  although 
the  deflections  are  relatively  close,  the  bending  moments  based  on  eq  15  do  not  approximate 
closely  the  actual  bending  moments,  especially  in  the  vicinity  of  the  load. 

Other  approximate  solutions  for  the  infinite  plate  were  discussed  by  Korunov”  in  1967.  As- 
suming that  Bernshtein's'  eq  10a  is  the  correct  expression  for  predicting  the  bearing  capacity, 
Korunov  proposed  the  empirical  expression  (for  h in  centimeters). 


r,  6.0 

P - — ah  in  tons 
cr  IOO 

or  rewritten 

P(r  00 ah2  in  kilograms  (22) 

and  then  showed  that  for  special  situations,  it  agrees  with  the  results  of  eq  10a.  Noting  that  eq 
22  is  based  on  of  24  kg/cm’,  it  follows  that 


* To  determine  the  stresses  under  the  load,  the  correction  derived  by  Westergaard  (ref.  142,  p.  275)  may  be 
used. 
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Figure  5.  A comparison  of  approximations  for  P 
with  exact  graph. 


I 


Figure  6.  Semi-infinite  plate  with  free  edge 
subjected  to  load  P. 


Note,  that,  according  toeq  22,  for  a given  u{,  PCf  is  proportional  to  the  second  term  in  eq  12a 
or  12b.  h \ whereas  eq  21,  derived  for  a 0,  is  proportional  to  the  first  term,  /r.  Note  also  the 
difference  between  eq  7,  suggested  by  Korunov,  and  eq  21,  derived  by  Panfilov.  A comparison  of 
various  approximate  expressions  for  Pcr  with  the  one  based  on  eq  10a  is  shown  in  Figure  5. 

It  appears  that,  instead  of  deriving  numerous  approximate  expressions  for  eq  1°  that  differ 
substantially  from  each  other  and  are  not  much  simpler  than  the  exact  expression,*  first  it  must  be 
established  whether  eq  10  is  suitable  for  predicting  the  bearing  capacity  of  floating  ice  plates  for 
loads  of  short  duration.  This  and  related  questions  will  be  discussed  later. 

Solutions  for  the  floating  semi-.nfinite  plate  with  a free  edge  subjected  to  latere'  toads  were 
presented  by  Westergaard  ’ in  1923,  Shapiro1”  in  1943,  and  Golushkevich"  in  1944,  using  Fourier 
integral  methods.  Shapiro's  results  were  verified  and  extended  by  Never1  in  1965. 

In  1950,  Zylev,1”  using  the  criterion  in  eq  8 presented  calculations  of  the  bearing  capacity  of 
a floating  semi-infinite  ice  plate  subjected  along  its  free  edge  to  vertical  and  horizontal  loads.  How- 
ever, Zylev  s approximate  solution  of  eq  9 for  the  vertical  load,  recently  included  in  a number  of  pub- 
lications,1* “ is  incorrect,  as  shown  below. 

For  the  se  i-inf inite  plate  shown  in  Figure  6,  Zylev'”  assumed  an  approximate  solution  of 
the  form 


w(x.  y)  = (cosh(ctx)  + T sinh(ax)lf(y) 


(23) 


where 


I’  = 1 for  x < 0 

(24) 

T = -1  for  x ->  0. 

Substituting  eq  23  into  differential  eq  9 with  q s 0,  he  obtained  an  ordinary  differential  equation  of 
fourth  order  for  f(y).  To  determine  the  four  constants,  he  used  two  regularity  conditions  at  infinity 
and  the  conditions 


• A prospective  user  of  eq  10  does  not  have  to  be  familiar  with  Bessel  functions  if  he  utilizes  the  C(a)  vs  a 
graph  shown  in  Figure  3. 
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My(X,  0)  0 


(25) 


(26) 


Note  that,  for  the  chosen  deflection  surface  (eq  23),  dw.'Bx  is  discontinuous  along  the  y-axis, 
which  is  not  the  case  in  an  actual  plate.  Th  > quantity  d\/dx3  is  also  discontinuous  along  the 
y-axis,  this  implies  that  for  the  assumed  deflection  surface  there  exists  a line  load  along  the  y- 
axis.  This  is  in  contradiction  to  the  assumed  plate  load  shown  in  Figure  6.  Furthermore,  along 
the  fiee  edge,  where  the  largest  stresses  are  anticipated,  the  boundary  conditions  for  a free  edge 
are  not  satisfied.  Therefore,  the  validity  of  Zylev’s  solution  for  the  semi-infinite  plate,  even  for 
the  determination  of  an  approximate  Pcf,  is  questionable. 

According  to  Zylev’s1”  results,  the  largest  bending  moment  takes  place  at  the  point  x 0 and 
y 1.14  \D.  On  the  basis  of  this  analysis 


Pct  O.HAd  e"A)‘lof)r 


(27) 


where 


0.2  tHb 
ST) 


(28) 


According  to  Saapiro’s  results,  omax  takes  place  under  the  load.  Utilizing  criterion  in  eq  8,  the 
load  at  which  the  first  crack  occurs  becomes 


P(r  S{a)oth2 


(29) 


where  S(«)  for  i 0.36  is  given  in  Figure  7. 

In  1960.  Panfilov’"  compared  the  values  of  the  load  PCJ  for  the  infinite  plate  as  well  as  the 
semi-infinite  plate.  The  corresponding  grap’is  are  shown  in  Figure  7.  This  comparison  shows  that 
P(r  for  the  semi-infinite  plate,  according  to  Zylev’”  (dashed  line),  is  much  higher  than  Pcr  accord- 
ing to  Shapiro1"  and  Golushkcvich.”  In  0 • h £ < 0,5,  it  is  even  higher  than  the  Prr  of  the  infinite 
plate.  In  view  of  this  comparison  and  the  obvious  errors  contained  in  Zylev’s  solution,  it  is  sug- 
gested that  eq  27  should  not  lie  used  for  the  analysis  of  the  semi-infinite  plate  with  a free  edge. 
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Figure  7.  Comparison  of  analytical  results.”  Figure  8.  (PCT)jnf  p|/(Pcr)semi  inl  pl  vs  b fj(. 
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According  to  Panfilov,’*  it  follows  from  Figure  7 that  for  0,07  p/jl  1.0 

<f>cr^inf  p!  ~2>  ^^cr^semi-inf  pi 
A more  precise  relationship  is  shown  in  Figure  8. 

On  the  basis  of  the  graph  for  the  semi-infinite  plate  shown  in  Figure  7,  Panfilov”  proposed 
for  the  interval  0,07  b l < 1 the  following  approximate  expression 


Panfilov104  attempted  to  derive  an  approximate  expression  for  P for  the  problem  shown  in  Figure 
6,  assuming  that 


*(*.  y)  wQex p-A(*  + 
\ 2 


(31) 


However,  the  result  obtained,  similar  in  form  to  eq  21,  is  of  questionable  value.  The  objections 
raised  in  connection  with  eq  21  also  apply  here.  Note  that  the  deflection  surface  (eq  31)  does 
not  satisfy  the  differential  eq  9 or  the  boundary  conditions  along  the  free  edge,  where  the  stresses 
are  determined  for  comparison  with  criterion  in  eq  8. 

The  semi-infinite  plate  subjected  to  equidistant  loads  P along  the  free  edge  was  analyzed  by 
Westergaard1’1  in  1923.  Similar  problems  were  solved  by  Panfilov101  ,0'  in  1963.  The  publications 
of  Shekbter  and  Vinokurova,1"  and  Korenev  and  Chernigovskaia”  also  contain  solutions  to  re- 
lated problems. 

The  solution  for  the  semi-infinite  plate,  simply  supported  along  the  straight  edge  and  sub- 
jected at  any  point  of  the  plate  to  a concentrated  force  P,  as  shown  in  Figure  9,  was  derived  by 
Kerr  in  19.)9.  Using  thi  method  of  images,  the  following  exact  closed  form  solution  was  obtained' 


y 


Figure  9.  Semi-infinite  floating  plate,  simply  sup- 
ported along  the  straight  edge  and  subjected  to  a 
load  P. 
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In  1971,  Palmer’6  utilized  this  solution  to  construct  a number  ot'  influence  surfaces  for  bend- 
ing moments. 

A numerical  solution  for  the  semi-infinite  plate,  f lumped  along  the  edge  and  subjected  to  a 
force  P at  a point  on  the  plate,  was  presented  by  Korenev6’  in  1960. 

An  analysis  of  a flouting  infinite  strip,  free  along  both  edges  and  subjected  to  a lateral  load, 
was  presented  by  Shapiro"6  in  1912,  utilizing  the  Fourier  integral  method.  Detailed  results  for 
similar  problems  were  presented  by  Panfilov16’  l0’  in  1966  and  1970. 

The  solution  for  a floating  infinite  strip,  simply  supported  along  both  edges  and  subjected  to 
a concentrated  force  P at  any  point  on  the  plate,  was  presented  by  Kerr’6  in  1959,  utilizing  the 
method  of  images.  The  resulting  deflection  was  given  as  a rapidly  converging  infinite  series  of 
fundamental  solutions  for  the  infinite  plate.  Other  solutions  for  this  problem  were  presented  by 
Westergaard"1  in  1923,  in  terms  of  Fourier  series  and  by  Nevel”  in  1965  in  terms  of  a Fourier  in- 
tegral. A solution  for  a similar  problem  was  presented  by  Panfilov16’  in  1966,  also  using  the 
Fourier  integial  method. 

The  infinite  strip,  with  clamped  boundaries,  was  analyzed  by  Nevel”  in  1965  and  by  Panfilov'  * 
in  1966,  using  Fourier  integral  methods. 

In  I960,  Kashtelian4*  presented  calculations  fc.  the  direct  determination  of  Pf  [that  is,  by 
eliminating  step  3 (p.  4)  in  the  above  procedure!  that  are  based  on  the  observation  that  the  carry- 
ing capacity  is  reached  when  the  wedges,  which  form  initially,  break  off.  However,  Kashtelian  s 
solution  for  the  wedge-shaped  plate,  on  which  his  calculations  are  based,  is  incorrect,  as  shown 
in  the  following. 


P / i 

) 


Figure  10.  A floating  rectangular 
comer  plate  with  free  edges  sub- 
jected to  load  P at  the  comer. 


Then,  utilizing  the  Bubnov-Galerkin  method,  for  a one-term  approximation  he  used 

6 o 

and  determined  from  it 


For  the  rectangular  comer  plate  with  free  edges,  shown 
in  Figure  10,  Kashtelian  assumed  an  approximate  solution 

w(x,  y)  - f exp!  -a(x  . y)lcos(ux)cos(ay)  (33) 

where  a and  f are  unknown  parameters.  From  the  condition 


oc  oc 

p \ l yw  dxdy 

"o  0 

Kashtelian  obtained 
f 4°2p 


(34) 


(35) 


(37) 
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Thus,  according  to  eq  35 

f JL.. 


It  should  be  noted  that  the  above  analysis  contains  an  error,  namely,  because  the  assumed  de- 
flection expression  (eq  33)  does  not  satisfy  the  lioundary  conditions  of  zero  moments  and  zero  she; . 
inn  forces  along  the  free  boundary,  eq  36  is  not  complete.  According  to  the  principle  of  virtual  dis- 
placements, the  proper  Bubnov-Galerkin  equation  for  a one-term  approximation  v/=  fWj(  x,  y)  is 


I I (D^w  v kwjw^dxdy 
J0  0 


[ M (x,  0)wt  (x,  0)dx 
Jo  'y 


& 


Fy(x,  0)w((x,  0)dx 


f 


♦ Mx(0,  y)w,  (O.y)dy  - f Vt(0,  y)w,(0,  y)dy  - Pwt((),  0)  = 0 

0 * 0 


(39) 


where  Mx  and  My 

iyo,y) 

Vylx.  0) 


are  Riven  in  eq  19  and 


- -D|w, 
-D|w, 
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(40) 


Comparing  the  / value  given  in  eq  38  with  the  corresponding  values  of  the  exact  solution  of  an 
infinite  plate,  / = P/(8vyD),  and  the  (incorrect)  approximate  solution  by  Zylev1”  for  a semi-infinite 
plate,  t = P (2\yD),  Kashtelian,4*  without  justification,  generalized  his  solution  for  the  rectangular 
comer  plate  to  a solution  for  a wedge  of  any  opening  angle  <f>  (Fig.  11)  by  assuming  that 


P 


Figure  11.  Floating  wedge  shaped 
plate  of  opening  angle  4>  subjected 
to  load  P at  the  tip. 


(41) 


an  equation  which  satisfies  eq  38  for  <t>  = e/2  and  the 
other  two  cases  (<f>  = 2n  and  <i>  * n)  mentioned  above. 
Utilizing  criterion  in  eq  8,  he  then  obtained  for  the  “fail- 
ure load”  of  a floating  wedge  plate  of  opening  angle  6 the 
expression 


(42) 


Note  that,  according  to  field  observations, 41  when  • 120°,  Pcr  - Pf , Thus,  according  to  eq  42, 
for  a floating  wedge  with  4>  = rr/2,  as  shown  in  Figure  10,  the  breakthrough  load  is 


P = (if  111— 

r \2/  0.966 


0.259  othZ. 


Observations  in  the  field  indicate  that  the  failure  mechanism  of  a semi-infinite  plate  subjected 
to  a force  P at  the  free  edfs  Proceeds  as  follows.  First,  a radial  crack  forms,  which  starts  under  the 
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Figure  12.  Failure  mechanism  nf  a Figure  13.  Failure  mechanism  for  a larpe  floating 

floating  semi-infinite  plate  subjected  plate  subjected  to  a load  P. 

along  the  free  edge  to  a load  P. 


load  and  propagates  normal  to  the  free  boundary.  This  is  followed  by  the  formation  of  a eircumfer- 
ential  crack  that  causes  final  failure,  as  shown  in  Figure  12.  According  to  Kashtelian, <a  the  failure 
load  for  this  case  is  equal  to  the  failuie  load  of  two  free  floating  wedges,  each  of  opening  angle 
<t>  n 2 


Pf  2 « 0.259 ofh'-‘  0.  r»lH /iL*.  (43) 

In  a similar  > ay.  Kashtelian4*  determined  the  Pf  for  an  infinite  plate.  Assuming  that  n is  the 
number  of  radial  cracks  and  that  the  n formed  wedges  are  all  of  equal  opening  angle,  i.t.,  <£n  = 2*  n, 
as  shown  in  Figure  13,  ttie  following  expression  for  the  failure  load  results: 


1 


n x 0.966 


o,h‘ 


Noting  that  n 2n  (j>n,  this  expression  may  also  be  written  as 


where  <t>  is  the  opening  angle  of  the  formed  wedges.  Note  that , with  decreasing  d>n.  the  load  Pf 
in  eq  44decreases  and  that  the  above  approach  does  not  take  into  consideration  the  effect  of  the 
wedge-in  moments  along  the  cracks. 

Kashtelian  showed  that  the  results  of  150  tests  on  floating  ice  plates  agree  closely  with  the 
bearing  capacity  values  based  on  eq  43  and  44.  In  view  of  the  errors  discussed  above,  however, 
this  agreement  is  not  convincing. 

An  approximate  solution  for  the  quarter  plate  with  free  edges  loaded  at  the  apex  was  also  pre- 
sented by  Westergaard1”  in  1948. 

An  exact  close  form  solution  for  the  quarter  plate  simply  supported  along  the  edges  and  sub- 
jected at  any  point  of  the  plate  to  a concentrated  force  P was  presented  by  Kerr”  in  1959,  using 
the  method  of  images. 

The  response  of  a narrow  infinite  wedge  resting  on  a liquid  base,  as  a beam  problem,  is  des- 
crib'd by  an  ordinary  differential  equation  with  a variable  coefficient.  This  equation  was  solved 
by  Dieudonee'’1  in  1957  by  means  of  the  Laplace  method  of  integrat  on.  Nevel”  in  1958  solved  it 
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Figure  14.  Wedge  shaped  plate  subjected 
to  a load  P = qb. 


lls in*;  the  method  of  Frobemus.  Nevpl’s  solution 
consisted  of  a sum  of  four  infinite  series  which 
were  evaluated  and  are  presented  as  graphs  in 
ref.  90.  An  approximate  solution  for  large  values 
of  x.  was  presented  by  Het^nyi.41 

An  early  attempt  to  determine  the  carrying 
capacity  of  a floating  ice  plate,  utilizing  a float- 
ing wedge  solution,  was  described  by  Papkovich"’ 
in  1915.  In  this  analysis  it  was  assumed  that  the 
wedge  response  is  governed  by  a modified  bending 
theory  of  beams  (Fig.  14)  by  stating  the  base 
parameter  as 


and  the  flexural  rigidity  El  as 


El<x) 


— (b 

12(1-,-')  \ 


k\x) 


}(b  * 2xtg?) 


(45) 


(46) 


where  y is  the  specific  weight  of  the  liquid.  The  term  (1  - c2 ) was  apparently  included  to  get 
plate  action  for  the  wedge.  The  deflection  was  assumed  in  the  form 

w(x)  Ae’Alcos(Ax) 

where 


and  the  unknown  constant  A was  determined  by  minimizing  the  total  potential  energy.  Substituting 
the  determined 


2A2P 

4'*  * 

into  eq  47  yields  the  deflection 

912  p . 

w(x)  — I e*Alcos(Ax). 

The  bending  moment  is 

U(x)  = -El  = -El(x)  2 A A2  e'*1  cos(Ax) 
dx2 


(49) 


(50) 
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aii.l  the  stresses  in  ttie  upper  and  lower  f.liers  were  obtained  as 


«(x) 


± MU) 
W(x) 


til  - •>  -1- 

A\‘  e **  cos(>x). 

1 -c” 


from  the  condition  do  dx  0,  the  position  of  the  largest  stress,  x n (4A),  was  determined.  Sub- 
stituting this  value  into  the  above  equation,  it  follows  that 


a 


max 


0.15 


A A2  Eh. 


(51) 


Utilizing  the  failure  criterion  in  eq  8,  emai  = o{.  it  follows  from  eq  51,  using  eq  48  and  49,  that 

P'-LT^V-  (52, 

Noting  eq  48,  the  above  expression  for  the  failure  load  of  a wedge  ot  opening  angle  4>  may  also  be 
written  as 


(53) 


Pointing  out  that  an  ice  plate  breaks  up  under  the  weight  of  an  icebreaker  into  wedges  and  that  P. 
in  eq  53  is  of  the  form  " f 


P,  = i4j  h2*  A2  h5/4  (54) 

Papkovich  suggested  that  eq  54  be  utilized  for  the  determination  of  an  empirical  expression  for 
the  breakthrough  load  of  an  ice  plate  by  determining  the  parameters  A t and  A £ from  field  test  data. 

Although  eq  53  is  only  an  approximation  (for  example,  the  corresponding  bending  moment  at 
x = 0 is  4 0),  its  dependence  upon  h is  identical  with  that  of  expressions  in  eq  12a  and  12b  for  the 
infinite  plate  and  eq  30  for  the  semi-infinite  plate,  respectively.  Even  the  term  bv/p£  appears  in 
the  proper  place.  This  observation  will  be  of  importance  in  the  discussion  of  test  data  presented 
in  ref.  98. 

For  solutions  to  other  plate  problems,  whose  response  is  governed  by  differential  eq  9,  refer- 
ence is  made  to  the  books  by  Schleicher,u’  Shekhter  and  Vinokurova,'29  Korenev,”  *'  and  Korenev 
and  Chernigovskaia”;  to  the  survey  articles  by  Korenev”  **  and  Savel’ev'”;  and  to  the  literature 
on  the  analysis  of  highway  and  airport  pavements. 

When  a floating  ice  plate  seals  the  liquid  base,  in  addition  to  the  buoyancy  pressure  kw(x,  y), 
the  liquid  exerts  a uniform  pressure  p*  on  the  plate.  In  such  cases,  an  additional  condition  has  to 
be  imposed  on  the  solution  to  reflect  this  situation.  The  unknown  p*  is  determined  from  this 
condition. 

If  the  assumption  that  the  liquid  is  sealed  and  incompressible  is  justified,  then  this  addi- 
tional condition  is 
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where  the  integration  extends  over  the  domain  of  the  plate  A. 

F loafing  plates  subjected  to  condition  expressed  in  etj  55  were  analyzed  by  Kerr”  ” and 
Nevel.  1 Kerr  and  Becker”  solved  plate  problems  by  assuming  that  the  sealed  liquid  is  compres- 
sible. They  showed  that  the  effect  of  the  sealed  liquid  depends  not  only  upon  its  relative  com- 
pressibility but  also  upon  the  sealed  volume  the  larger  the  sealed  volume,  the  smaller  the  seal- 
ability  effect.  This  result  suggests  that  the  use  of  eq  55  for  the  analysis  of  an  ice  plate  that 
covers  a river  or  a lake,  as  suggested  recently  hy  Mahrenholtz,*’  is  not  justified. 

The  analyses  reviewed  in  this  section  are  based  on  eq  9,  the  differential  equation  for  a homo- 
geneous and  isotropic  thin  elastic  plate.  In  an  actual  floating  ice  plate,  the  material  parameters 
vary  across  the  thickness  of  the  plate,  hence,  the  floating  ice  plate  is  nonhomogeneous.  This 
variation  is  very  pronounced  in  sea  ice  plates  as  well  as  in  a plate  whose  upper  surface  is  sub- 
jected to  very  low  air  temperatures. 

An  early  attempt  to  take  into  consideration  the  variation  of  Young’s  modulus  E (ref.  11,  p.  73) 
is  incorrect  because  the  investigators  did  not  take  into  consideration  that  when  E varies  across 
the  plate  thickness  the  resulting  stress  distribution  is  not  linear. 

According  to  recent  analyses  by  Newman  and  Forray,”  Assur,*  and  Panfilov,1"  when  Young’s 
modulus  E vanes  with  the  plate  thickness  h,  and  Poisson’s  ratio  i is  assumed  to  be  constant,  eq  9 
is  still  valid  if  the  flexural  rigidity  is 


0 


h-  z0 

z 2 E(z)dz 

'0 


and  the  position  of  the  reference  plane  is  determined  from  the  condition 


h - z0 


(56) 


(57) 


For  the  utilization  of  the  available  solutions  of  eq  9 also  for  nonhomogeneous  plates  with 
E = E(z),  it  had  to  be  shown  that,  except  for  eq  56,  the  corresponding  boundary  conditions  are  the 
same  as  those  for  homogeneous  plates.  This  was  done  recently  by  Kerr  and  Palmer,”  who  sys- 
tematically formulated  this  prohlem  utilizing  Hamilton's  principle  in  conjunction  with  the  three 
dimensional  theory  of  elasticity.  Kerr  and  Palmer”  also  showed  that  even  though  the  plane  sec- 
tion hypothesis  is  assumed,  the  resulting  bending  stress  distributions  are  not  linear  across  the 
plate  thickness.  An  example  is  given  in  Figure  15.  This  finding  suggests  that  the  well  known 
stress  equation 
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Figure  IS.  Stress  distribution 
in  the  plate  for  a given  F.(z). 
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Figure  16.  Deflection  vs  time 
curves  lor  a floating  ice  plate  and 
fixed  loads. 


utilized  by  various  investigators  in  conjunction  with  the  criterion  in  eq  8,  or  for  the  determination 
of  the  failure  stress  <>f  from  tests  on  floating  ice  beams,  may  not  be  applicable  in  general. 

It  should  be  noted  that  highly  concentrated  loads  often  cause  punch-tl  rough  failures,  for  such 
situations,  that  is,  for  very  small  a,  the  above  methods  may  not  be  suitable  and  a different  approach 
may  have  to  be  utilized  to  determine  the  breakthrough  load. 


METHOD  BASED  ON  VISCOELASTIC  THEORIES 

It  was  observed  in  the  field  that  for  loads  that  do  not  cause  an  instantaneous  breakthrough 
the  ice  plate  deforms  at  first  elastically  and  then,  with  progressing  time,  continues  to  deform  in 
creep,  especially  in  the  vicinity  of  the  load.  Two  characteristic  deflect  ion- vs-time  curves  for 
fixed  loads  P are  shown  in  Figure  16.  Curve  1 represents  the  case  when,  after  a time,  the  rates 
of  deformation  diminish  and  the  ice  plate  and  load  come  to  a standstill.  This  curve  corresponds  to 
a safe  load  for  any  length  of  time  under  consideration.  Curve  11  represents  the  case  when,  after  a 
time,  the  rates  of  deformation  increase  and  at  time  tf  the  load  breaks  through.  Thus,  the  load  that 
corresponds  to  curve  11  is  safe  for  time  t ■ ff,  but  then  it  has  to  be  moved  to  another  location  to 
prevent  breakthrough.  The  above  field  observations  suggest  that  for  an  analytical  determination 
of  breakthrough  loads  which  do  not  cause  immediate  failure  a viscoelastic  analysis  must  lie 
conducted. 

It  appears  that  the  small  deformation  theory  of  plates  may  be  sufficient  for  plates  which  fol- 
low curve  I.  However,  the  analysis  of  plates  which  respond  according  to  curve  11  is  more  compli- 
cated because  in  the  vicinity  of  the  load,  a region  of  prime  interest,  the  small  deflection  theory 
may  not  tie  valid  for  t approaching  tf.  Also,  as  the  plate  deflections  increase,  the  plate  may  start 
to  crack  - a phenomenon  not  predicted  by  the  usual  theories  of  viscoelastic  continua.  To  predict 
cracking,  a separate  failure  or  crack  criterion  muit  be  used.  Also,  after  the  first  crack  takes 
place  the  analysis  gets  even  more  involved  because  of  the  introduction  of  additional,  often  irreg- 
ular, plate  boundaries. 

For  an  analytical  determination  of  a "safe”  load  P < Pf  and  a "time  to  failure”  tf,  it  is  de- 
sirable to  have  one  viscoelastic  theory  for  floating  ice  plates  which  for  time  t - 0 yield  the  elastic 
response  and  for  t ' 0 yield  responses  according  to  curve  I or  11,  depending  upon  the  load  and  the 
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material  parameter5  of  the  ice  (which  in  turn  depend  upon  the  temperature,  distribution,  salinity, 
etc.).  This  theory  should  be  supplemented  by  a crack  or  failure  criterion  valid  for  the  elastic  and 

IdSt  1C  riiniiP,.  Thp  PlUQfir*  fhortfU  in  nnni<m»t;A..  _ . 


viscoelastic  range.  The  elastic  theory  in  conjunction  with  the  crack  criterion  a 
above  could  if  proven  correct,  he  a special  case  of  such  a general  theory.  mai 


* Oi 


discussed 


Anotoer  failure  criterion  was  pr .nosed  by  Zubov1"  in  19**2  and  by  Kobeko  et  «1."  ,n  1946  On 
the  basis  of  their  test  data,  they  concluded  that  for  loads  of  short  or  long  duration,  a floating  ice 
Plate  fails  under  the  load  when  a certain  deflection  wf  is  reached,  that  is  when 


w w. 

mai  f- 


(58) 


/ 


According  to  Kobeko  et  al.”  for  this  criterion  it  does  not  mat- 
ter whether  the  plate  deflections  are  purely  elastic  or  visco- 
plastic. as  shown  in  Figure  17.  The  criterion  in  eq  58  was 
also  adopted  by  Savel’ev  (ref.  121.  p.  438)  in  1963  for  the 
study  of  the  effect  of  temperature  and  salinity  on  the  carrying 
capacity  of  a floating  ico  cover. 

In  1961,  Panfilov”  proposed  the  above  criterion  for  float- 
ing ice  plates  that  are  cracked  in  the  dished  area.  His  justifi- 
cation was  that  in  the  dished  area  water  begins  to  flood  the 
upper  surface  of  the  plate,  with  a resulting  loss  of  base  pres- 
sure in  this  area.  It  may  be  added  that  the  flooding  of  the 
upper  surface  near  the  load  also  raises  the  temperature  of  the 
upper  layers  of  the  plate  to  about  O C,  thus  decreasing  the 
strength  of  the  ice  in  the  area  of  high  stresses. 

~'h  pu,e  6 ™ 
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Figure  17.  Illustration  ol  the 
failure  criterion  based  on  plate 
deflections. 
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where  w,  and  ft  are  given  in  centimeters.  In  this  connection  it  is  of  interest  to  note  that,  using 

criterion  wina  wf  ,n  conjunction  with  eq  59  and  the  solution  for  an  infinite  (uncracked)  elastic- 
plate  subjected  to  a concentrated  load  P 


B\7TT 


it  follows  that 


8y  yD 


= 2. 2V7T 


or 


P,  = 17.6 


120-e2)  ! 


ft2. 


160) 


Thus  according  to  the  criterion  in  eq  59.  the  breakthrough  load  Pf  is  proportional  to  ft2.  It  may  also 
be  of  interest  10  note  that  if  the  k'gest  deflection  of  the  plate  under  consideration  is  expressed  by 


the  equation 
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\Uieiv  i is  a coefficient,  then  a Pj  expression  of  the  form  shown  in  eq  54  corresponds  to  the 
criterion 


wt  * 7= r (61) 

\ n 

where  a and  are  coefficients. 

Test  data  are  needed  to  establish  whether  the  failure  criterion  in  eq  58  and  its  special  forms 
in  eq  59  or  61  are  indeed  valid  for  elastic  as  well  as  viscoelastic  deformations. 

In  the  early  attempts  to  take  time  effects  into  consideration  for  floating  ice  plates,  one  ap- 
prowih  utilized  the  solutions  for  elastic  bending  and  tried  to  fit  the  experimental  data  by  modify- 
ing the  elastic  constants  (ref.  11,  p.  53).  In  another  approach,  the  elastic  results  were  multiplied 
by  a time  factor  (1  , ut P),  where  t is  time  and  a and  /I  are  constants  to  be  determined  from  experi- 
mental data  (ref.  6,  eq  177).  However,  these  approaches  have  no  rational  foundation  and  their  re- 
sults are  of  questionable  value. 

Another  early  approach  was  based  on  Zubov's  hypothesis,  which  states  that  deflections  of 
ice  plates,  especially  at  comparatively  hi^h  temperatures,  are  caused  mainly  by  vertical  shearing 
forces  (ref.  154,  p.  49).  To  verify  Zubov’s  assumption,  Zvolinskii"6  analyzed  a plate  resting  on 
a liquid,  assuming  that  the  deformations  are  entirely  due  to  shearing  action  and  that  for  creep  de- 
formations the  material  obeys  Newton's  law  of  viscosity.  Although  the  resulting  differential  equa- 
tion was  relatively  simple,  because  of  ihe  prescribed  initial  conditions  the  obtained  solution  was 
rather  involved  Zvolinskii  (ref.  156,  p.  21)  stated:  "In  this  formula  the  result  is  nc:  self  evident, 
and  analyzing  it  does  not  help  us  to  visualize  the  picture  of  the  phenomenon.” 

Zvolinskii  used,  for  the  initial  condition,  the  elastic  detlection  surface  caused  by  shear  only. 
However,  according  to  some  experiments,  shortly  after  the  load  is  placed  the  deflection  surface 
agrees  closely  with  the  elastic  deflection  surface  due  to  binding  (ref.  8,  Fig.  18).  Also,  since 
the  elastic  deflections  are  relatively  small,  the  effect  of  assuming  that  the  elastic  deformations 
are  zero  seems  to  be  negligible  compared  with  the  introduced  error  of  assuming  shear  as  the  only 
force  responsible  for  creep  deformations.  This  assumption  was  made  by  Kerr, 49  who  attempted  to 
simplify  Zvolinskii's  analysis  in  order  to  study  the  characteristic  features  of  the  creep  deforma- 
tions based  on  Zubov's1'4  hypothesis. 

Recorded  observations  of  the  effect  of  static  loads  on  the  deformation  of  floating  ice  fields 
showed  (Fig.  16)  that  in  some  cases  the  rates  of  deflection  decreased  after  the  load  was  placed 
and  after  a certain  time  interval  the  plate  came  to  a standstill  (ref.  8,  p.  48,  ref.  154,  p.  146), 
whereas  in  other  cases  the  rates  of  deflection  increased  until  the  plate  collapsed  under  the  load. 
The  observed  decreasing  and  increasing  rates  of  deflection  should  result  from  a general  formula- 
tion of  the  problem.  However,  because  of  the  simplifying  assumptions  made,  it  was  necessary" 
to  incorporate  it  by  setting  up  two  separate  formulations  for  the  decreasing  and  increasing  rates 
of  deformation.  Although  some  of  the  results  obtained  did  agree  with  deflection  expressions  given 
by  Zubov  (ref.  154,  p.  24,  ref.  155,  p.  148),  because  of  tht  various  assumptions  made,  the  result- 
ing analysis  is  not  conclusive  for  the  determination  of  breakthrough  loads. 

The  assumption  that  the  predominant  deformations  of  a floating  ice  plate  are  caused  by 
shearing  forces  was  also  made  by  Krylov’4  in  1948. 
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The  intense  development  of  the  linear  theory  of  viseoelastieity  after  World  War  II  also  affected 
the  formulation  of  ice  plate  problems.  In  19-14,  Golushkevich”  presented  an  analysis  assuming 
that  ice  behaves  elastically  for  volumetric  deformations  and  viseoelastieally  for  deviatorie  defor- 
mations. His  formulation  was  base-1  on  the  linear  bending  theory  of  plates,  linear  constitutive 
equations,  and  the  assumption  that  tie  material  parameters  do  not  vary  across  the  plate  thickness. 

The  equations  obtained  were  linear.  The  special  ease  of  an  incompressible  material  was  analyzed 
in  detail. 

A general  formulation  for  viscoelastic  plates,  based  on  the  linear  bending  theory  of  plates  and 
the  assumption  that  the  constitutive  equation  is  a linear  relation  of  differential  operators,  was  pre- 
sented by  Freudenthal”  in  1958.  The  utilization  of  this  equation  for  floating  iee  plates  was  dis- 
cussed bv  Kheishin’6  in  1964.  As  a special  case.  Kheistiin  analyzed  an  infinite  ice  plate  sub- 
jected to  a concentrated  load  P , assuming  that  the  iee  is  incompressible  for  voluinetrie  deformations 
and  that  it  responds  like  a Maxwell  body  for  deviatorie  deformations.  A similar  problem,  when  the 
load  is  distributed  uniformly  over  a circular  area,  was  analyzed  in  1966  by  Nevel,"  who  also  pre- 
sented graphs  and  a comparison  w-'a  the  results  of  a test.  In  1970,  IAkunin”  presented  solutions 
for  various  load  distributions,  assuming  that  the  ice  responds  like  a four  element  model,  that  is,  a 
series  combination  of  a Maxwell  and  Kelvin  model.  In  the  above  analyses,  except  for  the  paper  by 
IAkunin,  it  was  assumed  that  the  material  parameters  are  constant  throughout  the  plate. 

As  discussed  before,  in  an  actual  floating  iee  cover  the  material  parametirs  vary  with  depth. 

In  an  attempt  to  take  this  into  consideration,  IAkunin4’  derived  an  approximate  formulation  for  a 
varying  modulus  of  elasticity  and  coefficient  of  viscosity,  and  solved  the  formulation  for  a variety 
of  load  distributions.  He  found  that,  as  in  the  elastic  case,  the  variation  of  material  parameters 
across  the  plate  thiekness  has  a profound  effect  upon  the  stresses  in  the  iee  cover. 

A viscoelastie  analysis  of  the  ice  cover  based  on  Reissner’s  theory  of  plates,  which  considers 
the  effect  of  bending  as  well  as  shearing  forces  upon  the  deformations,  was  presented  in  1967  and 
1968  by  Garbaccio.  Garbaeeio  assumed  that  ihe  ice  responds  like  a series  ' ombination  of  a Max- 
well and  Kelvin  model  and  that  the  material  parameters  are  constant  throughout  the  ice  plate. 

In  1961,  Panfilov,”  citing  shortcomings  of  linear  theories,  derived  a differential  equation  fur 
floating  ice  plates,  based  on  the  linear  bending  theory  of  plates  and  the  nonlinear  viscoelastic  con- 
stitutive equations  proposed  by  Voitkovskii.144  141  Additional  derivations,  along  the  same  line,  were 
presented  in  1970  by  Panfilov,'"  who,  however,  gave  no  solutions  to  the  derived  differential  equation. 

In  1962.  Cutliffe  et  al.,‘*  using  a nonlinear 
constitutive  equation,  made  an  attempt  to  analyze 
the  time-dependent  stresses  of  an  ice  eover. 

The  linear  bending  theory  and  a nonlinear 
constitutive  equation  were  also  used  by  Gar- 
baecio”  to  analyze  ice  plate  problems.  Gar- 
baccio attempteo  to  obtain  an  approximate  solu- 
tion of  the  resulting  nonlinear  formulation  by  a 
linearization  technique. 

In  the  absenee  of  reliable  analyses  for  pre- 
dicting the  bearing  capacity  of  iee  plates  sub- 
; 2 3 4 5 jeeted  to  loads  of  long  duration,  Panfilov,”  in 

t(  (hr)  6 1961,  constructed  from  field  test  data  the  graph 

Figure  18.  Breakthrough  loads  vs  breakthrough  s*lown  >n  P'Kure  18.  In  Figure  18,  t is  the  time 
time  (or  a floating  ice  cover  subjected  to  loads  Per>od  between  placement  of  the  load  and  break- 
o(  long  duration.  through,  Ff(0)  is  the  magnitude  of  the  load 
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just  sufficient  to  break  through  immediately  after  placement  on  the  plate  (at  tf  0),  as  discussed 
in  the  previous  section,*  and  P,(tf)  is  the  load  that  breaks  through  after  a time  t..  From  the  graph 
shown,  it  follows  that  Pj(tf)  < P^O)  for  tf  > 0.  Thus,  for  example,  a load  that  has  to  park  safely  on 
the  ice  plate  for  6 hours  should  be  smaller  than  0.4  Pf(0),  where  Pf(0)  is  determined  from  a separate 
analysis.  To  represent  analytically  the  graph  shown  in  Figure  18,  Panfilov  proposed  the  expression 

Pf(‘f)  1 

P,( 0)  1 f 0.75^  <62) 

where  t(  is  in  hours.  Solving  this  equation,  the  “safe”  storage  time  is  obtained  as  the  time  that 
is  smaller  than 


"Pf(°)  -Pf(rf)|3 
f [>75  Pf(0) 


(63) 


A graph  similar  to  the  one  shown  in  Figure  18  was  presented  and  discussed,  also  in  1961,  by 
Assur.1 

Korunov”  in  1968  pointed  out  that  eq  62  was  obtained  from  tests  on  ice  plates  under  specific 
conditions.  He  then  proposed  the  following  modification  of  eq  63 


L°- 


Pf(0)  - Pft S3 


75  Pf(0)n 


(64) 


where  K and  n are  correction  coefficients  which  take  into  consideration  the  shape  of  the  load  and 
the  outside  temperature. 

In  1970,  other  expressions  of  the  type  shown  in  eq  62  were  presented  and  discussed  by  Pan- 
filov.1" A related  discussion  is  presented  in  ref.  43. 


METHODS  BASED  ON  THE  YIELD  LINE  THEORY  OR  LIMIT  ANALYSIS 

The  yield  line  theory  was  utilized  for  the  analysis  of  continuously  supported  plates  by 
Johansen4’  in  1947  and  by  Bernell’  in  1952.  Persson1"  used  it  in  1948  for  the  analysis  of  a 
floating  ice  plate.  Assuming  that  the  yield  line  moment  per  unit  length  is  M0,  Persson  obtained 
for  the  case  shown  in  Figure  2 


Pt- 


(!+■/)(  1-0.62 a**) 


2/3>0‘ 


(65) 


Using  a similar  approach,  in  1961  Assur1  presented  for  the  breakthrough  load  the  expression 


(66) 


• In  the  previous  section,  it  is  denoted  for  brevity's  sake  as  P^. 
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The  method  of  Until  analysis  was  utilized  by  Meyerhof*’  in  I960  for  the  analysis  of  the  bear- 
ing capacity  of  float  mu  ice  plates.  Assuming  1)  that  the  ice  plate  is  thin,  rigid  and  ideally  plas- 
tic. 2)  that  it  can.  without  cracking,  resist  a full  plastic,  moment  M0.  and  3)  that  the  ice  obeys  the 
Iresca  yield  condition,  Meyerhof  obtained  for  the  case  shown  in  Figure  2 


0.0, r>  < a s 1.0. 


(67) 


Assuming  that  the  floating  ice  plate  before  failure  is  cracked  radially  into  numerous  wedges,  Meyer- 
hof obtained  for  the  same  case 


Pl 


4 n 


o 


M0 


0.2  < a < 1.0. 


(68) 


In  an  extensive  discussion  of  Meyerhofs  paper, ,J  Hopkins  questioned  the  degree  of  realism 
in  approximating  the  mechanical  behavior  of  ice  as  that  of  a rigid,  perfectly  plastic  material.  In 
the  same  discussion,  Wood  as  well  as  Hopkins  questioned  the  use  of  the  Tresca  yield  condition. 

Recently.  Coon  and  Mohaghegh1*  also  analyzed  the  floating  ice  plate  by  using  the  limit  analy- 
sis method  but  assumed  that  the  ice  obeys  Coulomb’s  law.  Fu  the  problem  shown  in  Figure  2 
they  obtained 


Pf  2rr(2.3  4 2.9a)M0. 


(69) 


Additional  results  and  discussions  were  presented  by  Coon  and  Mohaghegh1*  and  Meyerhof  •»  Re- 
lated results  were  published  by  Korenev*1  in  1955  and  Serebrianyi12*  in  1960. 


It  should  lie  noted  that  the  often  used  expression  for  the  limit  bending  moment  Un  = aJi2/4  is 
based  on  a stress  distribution  of  a homogeneous  plate,  as  shown  in  Figure  19a,  whereas  because 
of  the  thermal  gradient  in  the  plate,  the  distribution  of  limit  stresses,11*  assuming  that  a full  plas- 
tic moment  does  exist,  could  be  as  shown  in  Figure  19b,  Also,  the  assumption  that  the  ice  plate 

can,  without  cracking,  resist  a full  plastic  moment  UQ  may  not 
a be  realistic,  since  its  formation  was  not  observed  in  the  field. 

When  using  the  yield  line  theory,  it  may  be  more  realistic  to  work 
with  cracks  instead  of  yield  lines,  and  wedge-in  moments  instead 
of  the  plastic  moment  UQt  especially  along  the  radial  cracks 
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A comparison  of  the  various  Pf  expressions  presented  above 
with  PCf  formula  (eq  10a)  given  by  Bernshtein*  is  shown  in  Figure 
20.  For  comparison,  it  was  assumed  that  a,  = a and  that  Mn 
Oq/i  4.  Note  that  a different  number  in  the  demoninator  of  U 
shifts  only  vertically  a plotted  graph.  All  Pf/(oJ)2)  versus  a 
graphs  obtained  using  plasticity  methods  show  the  same  charac- 
teristics and  may  be  represented  by  a straight  line,  as  done  in 
eq  12a  or  12b. 


Figure  19.  Limit  stress  dis- 
tributions. a.  Homogeneous 
plate,  b.  Non  homogeneous 
piste. 
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Figure  20.  A comparison  ol  Pt  expressions 
with  PCf  lormul a (in  eq  10a). 


COMPARISON  OF  ANALYTICAL  AND  TEST  RESULTS 
General  Remarks 

The  mechanical  properties  of  ice  vary  drastically  in  the  vicinity  of  the  melting  (or  freezing) 
temperature  of  about  0CC.  Because  the  lower  surface  of  a floating  ice  plate  is  usually  at  the  melt- 
ing temperature,  the  plate  response  is  obviously  affected  by  it.  This  effect  is  especially  severe 
when  t he  upper  surface  is  also  subjected  to  near  0°C  temperatures,  because  then  the  temperature 
throughout  the  plate  is  approaching  the  melting  temperature. 

To  demonstrate  the  temperature  variations  with 
time,  consider  a floating  ice  plate  subjected  for  a long 
time  to  an  air  temperature  of  -10°C.  Assume  that  at 
time  t * 0 the  air  temperature  rises  to  -1°C,  the  cor- 
responding temperature  distributions  for  different  times 
are  shown  in  Figure  21.  Although  the  temperatures  at 
the  top  and  bottom  surfaces  are  constant  for  t > 0,  the 
temperatures  throughout  the  plate  vary  with  time.  Hence, 
if  two  identical  tests  are  performed  before  a thermal 
steady  state  is  established,  the  results  may  differ,  de- 
pending upon  the  time  (after  the  temperature  rise)  a par- 
ticular test  is  conducted. 

A similai  situation  takes  place  in  the  floating  test  beams  used  for  the  determination  of  the 
failure  stress  Oj,  because  after  a beam  is  cut  out  from  the  ice  the  side  walls  come  in  contact  with 
the  rising  water  and  the  outside  air. 

Another  thermal  problem  may  arise  in  a test  when  an  ice  cover  in  the  field  is  loaded  by  pump- 
ing water  into  a large  tank  that  rests  directly  on  the  cover,  for  then  the  bottom  of  the  tank,  which 
is  made  of  metal  or  canvas,  rests  on  the  ice,  and  the  upper  surface  of  the  ice  plate  in  the  contact 


Liquid 

Figure  21.  Temperature  distributions 
in  a Hoating  ice  plate  for  different 
times  0 < t < ». 
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region  is  subjected  to  near  0CC  temperatures.  This  type  of  loading  usually  causes  a change  in  the 
stress  distribution  and  a lowering  of  the  strength  of  the  ice  in  the  area  where  failure  usually  starts, 
thus  affecting  the  test  results, 

These  and  related  questions,  such  as  the  effect  of  a sharp  drop  of  the  air  temperature,  the 
rate  of  loading,  and  the  penetration  of  water  through  the  ice  plate  during  loading,  have  to  be  con- 
sidered when  the  test  data  of  floating  ice  plates  are  correlated  with  the  analytical  results.  In  the 
following,  various  test  results  are  discussed  and  correlated  with  analyses  presented  above. 

Effect  of  Bending  and  Shearing  Forces  on  Deflection  of  an  Ice  Cover 

As  shown  in  the  previous  sections,  an  analytical  determination  of  the  breakthrough  load  uti- 
lizes a formulation  for  the  ice  cover.  In  order  to  simplify  the  necessary  analyses,  such  a formula- 
tion contains  a number  of  assumptions.  It  is  essential  that  the  assumptions  made  be  justified, 
from  a physical  point  of  view,  since  otherwise  the  analytical  results  may  have  no  relevance  to  the 
actual  problem  under  consideration. 

One  such  assumption,  included  in  the  derivation  of  differential  eq  9,  states  that  a straight 
line,  normal  to  the  reference  plane,  remains  straight  and  normal  to  the  deformed  plane  (sometimes 
denoted  as  the  Kirchhoff  hypothesis).  Physically,  this  kinematic  assumption  implies  that  the  de- 
flections are  caused  hy  bending  stresses  only  and  that  the  effect  of  shearing  forces  is  negligible. 
This  assumption,  discussed  at  length  in  books  on  the  strength  of  materials,  has  been  proven  to  be 
justified  for  the  elastic  response  of  slender  beams  and  thin  plates  made  of  a variety  of  materials. 

On  the  other  hand,  basing  his  view  on  field  observations,  Zubov1”  in  1945  suggested  that  the 
deflections  of  an  ice  cover  are  mainly  ciused  by  shearing  forces,  and  hence  the  effect  of  bending 
upon  the  deflections  is  negligible. 

Because  the  resulting  equations  a-e  used  for  the  analytical  determination  of  Pf  (for  additional 
examples,  see  ref.  102),  it  is  essential  to  determine  whether  Kirchhoff’s  or  Zubov's  assumption  is 
to  lie  used  for  the  formulation  of  ice  cover  problems.  In  this  connection,  note  that  the  plate  de- 
flections due  to  a load  q,  which  is  distributed  over  a circular  area,  according  to  eq  9,  are1” 

w.(r)  = — TI  + ker,(a)ber(Ar)  - keiYa)beiIAr)l  0 < r < a 

y la  J = * 

(70) 

w2(r)  =.  — |ber  (a)ker(Ar)  - bei  (a)kei(A.r)j  a < r < «, 

where  A *=  \y/D,  whereas  the  differential  equation  for  an  ice  plate,  according  to  Zubov’s  hypothe- 
sis, is 

Gh  y2  w - yw  = - q (71) 

where  G is  the  shearing  modulus  and  the  corresponding  deflections  are 

Wj(r)  = [l  - (Ka)Kj(«8)/0(Kr)J  0<r<a 

w2(r)  = S (xa)/j(Ka)K'0(/cr)  a < r < «• 


(72) 
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Figure  22.  Comparison  of  plate  deflection  curves  based  on 
bending  and  shear  theories. 
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Figure  23.  Comparison  of  ice  plate  deflections 
due  to  loads  of  short  duration  at  -15°C  < T < -7°C.‘ 
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Figure  24.  Comparison  of  ice  plate  deflections  due  to 
loads  of  short  duration  at  0°C. 


where  lQ,  KQ,  / x and  fC  t are  Bessel  functions  and  * = \ y (Gh) . To  show  the  different  nature  of 
the  deflection  curves  based  on  these  two  assumptions,  eq  70  and  72  were  evaluated  numerically 
for  h = 10  cm,  v = 0.3,  and  E 10.500  kg/cm'.  For  a/ft  = 1 aad  5,  the  corresponding  value  of  G 
was  determined  using  the  condition  that  the  largest  deflections  w(0)  for  both  theories  are  equal. 
The  results  are  shown  in  Figure  22. 

Note  that  the  response  of  an  ice  cover  according  to  Zubov'”  is  identical  to  the  response  of 
the  shear  layer  in  the  Pasternak  foundation.” 
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As  early  as  1929,  Bernshtein*  compared  the  deflections  of  an  ice  field  on  the  Volga  River  sub- 
jected to  loads  of  short  duration,  at  air  temperatures  of  -1  VC  • T -7  C,  with  cot  res  ponding  re- 
sults based  on  eq  9 This  comparison  is  shown  in  Figure  23.  Since  the  agreement  is  very  close, 
it  was  concluded  that  the  use  of  eq  9,  and  hence  Kirchhoff’s  hypothesis,  is  justified  for  the  formu- 
lation of  ice  plate  problems  subjected  to  loads  of  short  duration. 

in  1968,  Shmatkov'"  compared  test  data  of  an  ice  plate  on  Lake  Baikal  subjected  to  a vertical 
load  of  short  duration  but  at  air  temperatures  of  about  0 C with  analytical  results  based  on  eq  9 
and  71.  This  comparison  is  given  in  Figure  24.  On  the  basis  of  these  data,  Shmatkov  concluded 
that  at  air  temperatures  of  about  0'  C the  deformations  are  mainly  caused  by  shearing  forces. 

This  conclusion  raises  a serious  question  about  the  effect  of  the  air  temperature  upon  the 
range  of  validity  of  ec  9 and  71  for  the  formulation  of  ice  covers.  A comparative  study  involving 
more  test  data,  especially  at  air  temperatures  near  0 C,  is  urgently  needed  to  clarify  this  important 
question.  In  these  tests,  a special  effort  should  be  made  to  separate  the  elastic  from  the  nonelastic 
deformations.  It  may  also  be  advisable  to  note  the  difference  between  the  crystallographic  structure 
of  an  ice  cover  formed  over  a lake  in  which  the  water  is  essentially  at  rest  and  that  over  a river  in 
which  the  water  moves  at  a certain  velocity,  and  the  effect  of  a different  crystallographic  structure 
upon  the  mechanical  properties  of  an  ice  cover. 


Deterninatioo  of  Pf( 0) 


Test  results  and  their  relationship  to  the  allowable  load  given  by  the  analogy  method  were 
discussed  by  Kliucharev  and  lziumov”  in  1943  and  bv  Kobekoet  al.”  in  1946.  In  1960,  Gold** 
compared  eq  4 with  the  field  results  of  the  Canadian  pulp  and  paper  industry.  The  conclusion 
from  this  comparison  was  that  the  formula  given  in  eq  4 is  not  sufficient  for  the  determination  of 
failure  loads  since  the  presence  of  cracks,  thermal  stresses  and  natural  variation  in  effective 
thickness  is  not  considered.  Another  reason  could  be  that  the  failure  load  Pf  is  not  proportional 
to  h 2 hut  may  be  a more  complicated  function  of  h,  as  indicated  by  eq  10  and  12.  Additional  re- 
sults were  presented  by  Gold”  in  1971. 
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Figure  25.  Laboratory  test  results  for  infinite 
plates” 


In  order  to  establish  which  of  the  various 
formulas  for  Pcr  and  P^O)  obtained  using  the 
criterion  cmax  = o(  are  suitable  for  predicting 
the  carrying  capacity  of  a floating  plate  sub- 
jected to  loads  of  short  duration,  in  the  follow- 
ing the  analytically  obtained  Pf(0)  values  are 
compared  with  corresponding  results  from  tests 
conducted  on  floating  ice  plates. 

Since  the  analyses  are  based  on  an  elastic 
theory,  only  the  results  of  tests  with  very  short 
loading  times  to  failure  are  of  interest.  Such 
tests  were  recently  conducted  by  Panfilov”  in 
the  laboratory  as  well  as  in  the  field.  The  lab- 
oratory tests  were  conducted  at  -1(TC.  The 
floating  plate  was  loaded  by  means  of  stamps 
of  the  dimensions  shown  in  Figure  25.  The 
load  for  each  test  was  placed  statically  at 
rates  which  caused  breakthrough  within  5 to 
20  sec.  In  addition  to  the  failure  loads  Pf, 
loads  at  which  the  first  radial  crack  occurred 
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PCT  were  recorded.  The  laboratory  tests  were  conducted  with  fresh  and  salt  water  ice.  The  ice 
plate  thickness  varied  from  7 to  30  mm.  The  field  tests  were  conducted  on  thicker  ice  plates.  The 
ice  plate  was  loaded  by  placing  metal  water  tanks  on  a structure  which  in  turn  rested  on  the  ice 
plate  and  simulated  the  contours  of  wheel  loads.  The  ice  strength  af  was  determined  from  floating 
cantilever  tests  with  the  load  acting  downwards.  Additional  details  are  contained  in  ref.  98. 

The  results  of  56  laboratory  tests  for  the  infinite  plate  are  shown  in  Figure  25. 

The  failures  followed  the  usual  pattern  first,  the  formation  of  radial  cracks  that  emanated 
from  the  region  under  the  load,  then,  the  formation  of  circumferential  cracks,  at  which  time  the  load 
broke  through  the  plate. 

In  Figure  25,  Curve  1 is  the  P,r  (o{h~),  according  to  the  analyses  by  Bernshtein,"  Colushkevich,” 
and  Wyman.1”  Curve  11  is  proposed  hy  Panfilov”  as  representing  the  test  data  and  is  described 
by  the  equation 

— 1.25  * 1.05  ^ . 

0{h~  l 

This  equation  was  ohtained  by  an  averaging  process.  The  test  data  show  a scatter  in  a relatively 
narrow  band. 

Before  proceeding  with  the  discussion  of  these 
test  results,  m the  following  a different  concept  for 
the  evaluation  of  ice  plate  tests  is  introduced.  This 
is  necessary  because  averaging  curves,  such  as 
curve  II,  are  not  suitable  for  most  engineering  pur- 
poses. 

From  an  engineering  point  of  view,  there  is  a 
need  to  determine  safe  loads,  at  which  an  object 
may  move  or  park  briefly  on  a floating  ice  plate, 
or  breakthrough  loads,  for  the  design  of  ice  break- 
ers, at  which  the  plate  definitely  collapses.  These 
loads  may  be  obtained  by  introducing  into  the  re- 
sults of  field  tests  an  upper  envelope  U and  a lower 
envelope  L,  as  shown  in  Figure  26.  It  is  reasonable 
to  expect  that  the  area  under  envelope  L contains 
safe  loads  and  the  area  above  envelope  U the 
breakthrough  loads.  The  area  between  the  envelopes  is  the  region  of  the  test  failure  loads  and 
nothing  definite  can  be  said  about  it  with  r^pect  to  safety  or  breakthrough.  Thus,  only  the  re- 
gions ahove  curve  U and  below  curve  L are  of  interest  and  the  test  results  are  needed  to  separate 
these  two  regions. 


Brcan  through  l OOJ^ 


P,  •» 


Soli-  t O'U 


Figure  26.  Illustration  of  new  concept 
for  evaluating  breakthrough  loads  from 
test  data. 


For  the  test  da'a  of  infinite  plates  showt  in  Figure  25,  the  upper  envelope  U may  be  repre- 
sented by  the  equat.on 


and  the  lower  envelope  L by  the  equation 
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1.0  ♦ 1.2-  • 
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Therefore,  if  the  hounds  shown  in  Figure  25  should  prove  reproducible  by  other  investigators,  a 
safe  load  (for  loads  of  short  duration  and  T -10  C)  could  be  determined  from  the  condition 


where  a{  is  obtained  from  a floating  cantilever  beam  test  loaded  downward. 

According  to  test  data  stiown  in  Figure  25 

(Pf'PSt)1  = 2 Pcr- 

Note,  however,  that  the  values  for  these  two  cases  are  usually  not  the  same, 

Panfilov  observed  that,  if  P(,f  is  ttie  load  at  which  the  first  crack  takes  place,  then 

ptest  ~ - piesi 
cr  3 rf 


(73) 


(74) 


(75) 


From  the  atrove  two  equations,  it  then  follows  that 


ptest  ~ * p 
cr  3 cr' 


(76) 


A proper  analysis  should  yield  that  PCf  is  equal  to  P**\  Possible  reasons  that  this  is  not  so  in 
eq  76  are’  l)  The  af  values  used  in  Figure  25  are  those  obtained  by  loading  the  cantilever  beam 
downward,  whereas  for  the  determination  of  PCf  the  tensile  stresses  that  crack  the  plate  are  in  the 
lower  fibers  of  the  plate  where  of  is  smaller  because  of  the  higher  temperatures'.  2)  The  stress 
distribution  is  not  linear  across  the  plate  thickness  and  the  stresses  in  the  upper  fibers  are  lar- 
ger than  those  in  the  bottom  fibers,  whereas  the  analyses  and  test  evaluation  are  based  on  a linear 
distribution  with  equal  stresses  at  the  top  and  bottom  fibers;  and  3)  The  criterion  a = a.  mav  nr 
be  valid. 

According  to  analytical  results  by  Kashtelian,*'  for  an  infinite  plate  that  cracks  into  five 
wedges  (<j!>n  = 2rr/5) 


P,  <t>  n 
——  = 2.08—  ~ 0.8 
oth2 


and  when  plates  crack  into  six  wedges  (<£n  = w/3) 


Ofh2 

Thus,  according  to  this  analysis.  Pf  values  are  obtained  which  are  far  below  the  test  data  pre- 
sented in  Figure  25. 


30 


THE  BEARING  CAPACITY  OF  FLOATING  ICE  PLATES 


0 8‘ 


p 0 «f; if' . •*  . 

V **«•  ♦>*  ViT 


02 


0 6 0 8 I 0 


4 V 

v /, 

Figure  27.  Results  lor  semi~in finite  plate.1 
(Symbols  as  in  Fig.  25 .) 


Also  compare  the  graphs  presented  in  Figure 
if)  with  the  test  data  of  Figure  25.  Note  that  the 
upper  graphs  in  Figure  20  are  based  on  UQ  oj i2/t 
and  that  they  may  be  shifted  toward  the  test  data 
hy  choosing  a larger  number  in  the  denominator  of 

M0- 

The  test  results  for  a seini-infimte  plate 
subjected  to  an  edge  load,  as  shown  in  Figure  6, 
are  presented  in  Figure  27.  The  failures  followed 
the  usual  pattern  first,  the  formation  of  a c rack, 
which  emanates  under  the  ioad  and  is  normal  to 
the  free  houndary;  tl-  u the  formation  of  a circum- 
ferential crack  at  which  the  two  wedges  break  off. 

In  Figure  27,  Curve  1 is  the  P , according  to 
the  analyses  of  Shapiro1”  and  Golnshkevich.” 

Curve  11  was  proposed  hy  Panfilov’'  as  represent- 
ing the  test  data,  which  show  a scatter  in  a rela- 
tively narrow  band.  It  is  described  by  the  equation 


P h 

-2—  , 0. 15  , 0.38  l 

■ g 1 

<if  n 


It  can  be  easily  verified  that  the  upper  envelope  U is  described  by  the  equation 


0.58  . 0.27  £ 

l 


and  tho  lower  envelope  L by  the  equation 


= 0.35  c 0.30  - . 

i 


Hence,  if  the  bounds  shown  in  Figure  27  should  prove  to  be  reproducible  by  other  investigators,  a 
safe  load  for  the  crossing  of  a long  gap  in  a floating  plate  (a  bridge  between  two  semi-infinite 
plates)  could  be  determined  from  the  condition 

P [0.  35  c 0.30 -A_\ 

V \'/Ty/ 

On  the  other  hand,  the  breakthrough  load  for  a semi-infinite  plate,  often  needed  for  the  design 
of  icebreakers,  should  satisfy  the  condition 

p ■ ' oa^r)‘’'"2 

whore  of  is  determined  from  a floating  cantilever  test  loaded  downward. 

According  to  tho  test  data  shown  in  Figure  27,  for  0.1  < fo/£  <,  1.0 
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10  ^ 

Panfilov  observed  that  also  for  the  semi-infinite  plate 
ntest  ^ 2 test 

Ct  3 *f  • 


(77) 


(78) 


f 10m  eq  77  and  78,  it  then  follows  that 


plest 

cr 


1.1  P 


cr' 


(79) 


In  view  of  the  tliree  possible  shortcomings  listed  in  the  discussion  of  the  infinite  plate  this 
agreement  is  very  close. 


Panfilov's  test  results  for  the  infinite  and  semi-infinite  plate  show  that 

(p/esl)  r 2.7  /p,est\ 

\ /inf  |>la(e  y f /semi-inf  plate" 


(80) 


This  contradicts  the  -xperimental  findings  reported  by  Kashtelian  (ref.  48.  p.  33).  Equation  80 
indicates  that  the  effect  of  the  wedge-in  moments  is  not  negligible  if  one  attempts  to  compute  P 
analytically  from  wedge  solutions.  Without  the  wedge-in  moments,  P(  of  the  infinite  plate  wouW 
be  equal  to  twice  the  Pf  of  the  semi-infinite  plate.  In  this  connect  ior,  nate  the  corresponding  re- 
lationship obtained  analytically  and  shown  in  Figure  8. 

According  to  Kashtelian,"  for  the  observed  wedge  formation,  for  a semi-infinite  plate 
<6  = n '2  and  hence 


Pcr/(offT)  0.518 


a value  which  agrees  with  the  test  data  shown  in  Figure  27  for  a < 0.4. 

Other  test  data  for  loads  of  short  duration  were  obtained  by  lAkunin,"  however,  these  results 
were  not  available  for  review. 


Determination  of  Pf(i) 

r « fW  iCC  C°VPrS  subjec,ed  t0  loads  of  long  duration  were  reported  in  refs. 

’ 8.  45,  58  and  59.  More  recent  test  results  were  presented  by  Sundberg-Falkenmark,'"  Franken- 
stein, Panfilov,  1 * 111  Stevens  and  Tizzard,'”  and  lAkunin.45  Although  some  writers  compared 
their  test  data  with  analytical  results  and  found  satisfactory  agreement  for  certain  situations,  a 
systematic  study  of  available  test  data,  supplemented  with  new  test  results,  is  needed  to  est’ab- 
lsh.  first  the  proper  plate  theory  for  ice  covers  which  will  predict  the  deflections  as  a function  of 
time,  and  then  a failure  criterion  for  the  determination  of  Pf(t)  and  tf. 

In  connection  with  the  above  studies  it  may  also  be  useful  to  note  the  test  results  presented 
by  Blacky  Brunk,'5  Butiagin."  Frankenstein,"  Gold  et  al.,“  Korzhavin  and  Butiagm,"  and 
Shishov,1*1  as  well  as  the  discussions  by  Assur,5  Dykins,”  lAkunin,"'  Pister,uriand  the  discus 
sion  in  ref.  140. 

Determination  of  a( 

For  the  analytical  determination  of  Pf(0),  the  value  af  is  needed.  It  is  determined  usually 
from  a beam  cut  out  from  an  ice  plate  and  tested  in  situ.  A detailed  description  of  such  tests 
was  given  by  Butiagin  (ref.  16,  section  IV).  A cantilever  test  beam  is  shown  in  Figure  28. 
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Figure  28.  Cantilever  text  beam  (or  the  determination  o(  u.. 


Other  test  data  were  presented  by  Brown,"  Frankenstein,”  " Sokolnikov,"*  Tabata  et  al.,‘" 
Tauriainen,1"  and  Weeks  and  Anderson.147  Related  questions  were  discussed  by  Butiagin," 
Frankenstein,"  Kerr  and  Palmer,”  Lavrov,”  71  Peschanskn,11*  Savel’ev,"1  Smirnov,"’  and  Weeks 
and  Assur.141' 

Ir  order  to  establish  a standard  procedure  for  the  determination  of  a(,  it  should  be  of  interest 
also  to  determine  the  effect  of  the  rate  of  loading  upon  o(,  as  well  as  to  clarify  why  Franken- 
stein,” ” using  the  test  setup  shown  in  Figure  28,  found  that  the  determined  of  value  was  higher 
when  P acted  upwards,  whereas  Butiagm,14  using  the  same  setup,  reported  that,  according  to  his 
test  results,  the  of  value  was  higher  when  P acted  downwards. 


SUMMARY  AND  RECOMMENDATIONS 

When  utilizing  a floating  ice  plate  for  storage  purposes  or  as  a pavement  for  moving  vehicles, 
there  is  a need  to  know  the  magnitude  of  the  breakthrough  load  Pf(t)  and  the  corresponding  time  to 
failure  tf.  Until  now,  there  has  been  no  general  theory  in  the  literature  suitable  for  the  prediction 
of  Pf(t).  The  majority  of  papers  on  the  bearing  capacity  of  ice  plates  have  dealt  with  the  deter- 
mination of  Pf(0),  that  is,  the  load  which  is  just  sufficient  to  break  through  the  ice  immediately 
after  it  is  placed  on  the  ice  cover.  Only  a few  papers  have  dealt  with  the  determination  of  Pf(f). 
The  procedures  for  the  determination  of  Pf(0)  and  Pf(t)  are  summarized  in  the  following: 


P/0) 

P/t) 

[ Based  on  elasticity  analyses 

Based  on  plasticity  analyses 

Analogy  method 
for  determina- 
tion of  Pall. 

Determination  of 
Pc  r based  on 
elasticity  theory 
of  plates  and  cri- 
terion <7mx=of. 
Then  correlation 
of  Pcr  and  Pf(0). 

Direct  determin- 
ation of  Pf(0)  by 
analyzing  the 
cracked  plate. 
Use  of  elasticity 
theory  and  cri- 
terion  aml=af. 

Determination 
of  Pf(0)  using 
yield  line 
theory. 

Determination  of 
P(0)  using  limit 
load  theory. 

Lse  of  a visco- 
elastic theory 
in  conjunction 
with  a failure 
criterion. 

Attempts  to  determine  P^O)  are  based  on  elasticity  as  well  as  placiicity  theories 
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The  foundations  of  the  ana  log  method,  wltieh  utilizes  relationships  of  an  elasticity  analysis, 
are  quest lonahle.  Thus,  the  results  obtained  with  this  method,  although  very  simple,  should  be 
used  with  caution.  In  this  connection,  note  the  position  of  eq  7 in  Figure  25  as  compared  wit ti 
some  findings  by  Gold.’4  ” 

Another  approach  is  based  on  the  elasticity  theory  ol  plates.  In  this  procedure,  for  the  given 
load  the  maximum  stiess  in  the  plate  <>max  is  determined  first,  then  eq  8,  uniax  ul,,  is  used  to  deter- 
mine Pcr,  a load  which  is  just  sufficient  to  cause  the  first  crack.  Since,  according  to  field  tests  for 
infinite  and  semi  infinite  plates,  Pf(0)  Pcr,  an  empirical  relation  between  PCTand  Pf(0)  is  needed  for 
the  determination  of  Pf(0).  Equation  74,  which  is  based  on  data  by  Panfilov!*'  if  proven  to  be  gen- 
erally valid,  could  be  used  as  such  an  empirical  relation  for  the  infinite  plate.  In  this  procedure,  o. 
is  determined  from  a floating  ice  beam  that  fails  in  tension  in  the  bottom  region  of  the  cross  section. 

In  still  another  approach,  the  empirical  relation  is  eliminated  and  Pf  (0)  is  determined  directlv, 
by  using  the  elasticity  theory  for  the  analysis  ol  the  cracked  ice  plate,  which  consists  of  wedges  ' 
that  emanate  from  the  loaded  region  and  the  assumption  that  Pf( 0)  is  reached  when  the  wedges  break 
off.  Equation  8 is  also  utilized  as  the  crack  criterion.  The  value  f f is  obtained  from  a floating 
ice  beam  that  fails  in  tension  in  the  upper  region  of  the  cross  section 

The  publications  that  follow  either  of  these  two  approaches  contain  several  questionable  as- 
sumptions, for  example,  although  in  a floating  ice  plate  tlie  material  parameters,  especially  E, 
usually  vary  throughout  the  thickness,  the  expression  valid  for  a linear  distribution  of  bending 
stresses  is  used  exclusively  for  the  determination  of  the  maximum  stress  in  the  plate.  Also,  the 
use  of  above  equations  for  the  determination  of  a,  from  a beam  test  may  not  be  justified.” 

Another  questionable  practice  is  the  utilization  of  eq  8 as  the  failure  criterion.  Equation  8 is 
the  well  known  maximum  stress  criterion.”  ,M  It  implies  that  the  failure  stress  af  is  not  affected 
by  any  other  stresses  at  the  point  of  failure.  Tests  have  shown  that  eq  8 is  applicable  to  a variety 
of  brittle  materials  when  not  subjected  to  hydrostatic  compression.  Although  many  publications 
dealing  with  the  bearing  capacity  of  floating  ice  plates  utilize  eq  8,  not  a single  publication  could 
be  located  which  describes  test  results  that  prove,  or  disprove,  the  validity  of  this  criterion  Tor 
floating  ice  plates.  This  situation  is  very  unsatisfactory,  since  omax  in  plates  is  usually  biaxial, 
whereas  '.he  of  value  is  determined  from  a test  with  uniaxial  bending  stress.  In  1970,  Panfilov110 
suggested  the  criterion 

ol-,u,2<Of  (81) 

which  is  the  two-dimensional  version  of  the  well  known  maximum  strain  criterion.  However,  Pan- 
filov110 did  not  offer  sufficient  experimental  data  to  justify  the  use  of  this  criterion  either.  In  the 
literature  on  the  mechanics  of  materials,  several  other  failure  criteria  are  described  that  may  or  may 
not  be  suitable  for  floating  ice  plates.  For  an  early  discussion  related  to  plates  on  a Winkler  base, 
the  reader  is  referred  to  section  9 of  ref.  123.  It  appears  that  first  it  has  to  be  established  whether 
the  simple  criterion  in  eq  8,  which  is  also  applicable  for  materials  with  different  cf  values  for  ten- 
sion and  compression,  is  valid  for  floating  ice  plates  subjected  to  vertical  and  in-plane  loads. 

An  additional  shortcoming  of  the  publications  that  analyze  the  cracked  plate  is  that  the  investi- 
gators neglect  the  wedge-in  moments  in  the  radial  cracks.  This  does  not  seem  to  be  permissible, 
in  view  of  the  tests  by  Panfilov,”  who  found  that  Pf(0)  of  an  infinite  plate  is  larger  than  2Pf(0) 
for  a semi-infinite  plate. 

The  approaches  for  the  determination  of  Pf(0)  that  are  based  on  plasticity  theories  utilize  the 
yield  line  or  limit  load  analysis.  For  a discussion  of  a possible  shortcoming  of  these  two  analyses, 
the  reader  is  referred  to  the  listed  references.  Note  that  the  yield  line  theory  is  conceptually  related 
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to  the  approach  discussed  above,  which  analyzes  the  cracked  plate.  In  this  connection,  it  may  be 
more  realistic  to  work  with  cracks  instead  of  yield  lines  and  wedge-in  moments  instead  of  plastic 
moments,  especially  along  the  radial  cracks. 

In  view  of  the  variations  of  ice  properties  in  an  actual  ice  cover  and  their  effect  upon  Pf(0), 
it  may  be  advisable  from  a practical  point  of  view  to  use  the  concept  presented  in  Figure  26.  Its 
theoretical  justification  is  that  the  straight-line  upper  or  lower  bounds  of  Pf(0)  are  of  the  form 

P, 

= A t Ba 

of  h 2 

which  relates  the  concept  to  the  various  analyses  discussed  above.  This  approach,  if  restricted  to 
straignt-line  bounds,  is  essentially  the  same  as  the  one  discussed  by  Papkovich,11’  except  for  the 
introduction  of  the  notion  of  upper  and  lower  bounds  for  PftO).  Also,  note  the  similarity  of  the  trend 
of  the  graphs  and  test  data  shown  in  Figures  4,  20  and  25. 

The  experimental  data  for  Pf(0)  presented  by  Panfilov*'  (Fig.  25  and  27)  show  little  scatter. 
More  test  data  are  needed  to  establish  whether  the  Pf  values  for  other  ice  plates,  tested  under  dif- 
ferent conditions,  fall  in  the  same  range. 

The  analytical  determination  of  Pf(t)  has  received  much  less  attention  than  the  determination 
of  Pf(0).  It  is  reasonable  to  assume  that  ttie  necessary  formulation  consists  of  a viscoelastic  plate 
theory  and  a failure  criterion.  Thus,  it  is  essential  first  to  establish  the  range  of  validity  of  a 
simple  formulation  consisting  of  a linear  viscoelastic  plate  theory  (a  bending  theory,  a shear  theory, 
or  a combination  of  both  effects)  in  conjunction  with  a failure  criterion  of  the  type  expressed  in 
eq  58. 

Until  reliable  analytical  methods  are  developed  for  predicting  Pf(t)  and  tf,  from  a practical 
point  of  view,  it  appears  advisable  to  establish  whether  the  empirical  relation  expressed  in  eq62 
or  a similar  expression,  as  proposed  by  Assur,2  is  generally  valid.  The  test  results  needed  for  this 
purpose  are  also  necessary  for  formulating  the  proper  failure  criterion  as  well  as  for  establishing 
the  validity  of  a choseii  viscoelastic  plate  theory. 
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